Shear-induced structuration of confined carbon black gels: 
Steady-state features of vorticity-aligned fiocs 
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Various dispersions of attractive particles are known to aggregate into patterns of vorticity-aligned 
stripes wlien slieared in confined geometries. We report a thorough experimental investigation of 
such shear-induced vorticity alignment through direct visualization of carbon black gels in both 
simple plane shear and rotational shear cells. Control parameters such as the gap width, the strain 
rate, and the gel concentration are systematically varied. It is shown that in steady states the 
wavelength of the striped pattern depends linearly on the gap width h while being insensitive to 
both the gel concentration C and the shear rate 7. The width of the vorticity-aligned floes coincides 
with the gap width and is also independent of C and 7, which hints to a simple picture in terms 
of compressible cylindrical floes. Finally, we show that there exists a critical shear rate 7c above 
which structuration does not occur and that 7c scales as h~" with a = 1.4 ± 0.1 independently of 
C. This extensive data set should open the way to quantitative modelling of the vorticity alignment 
phenomenon in attractive colloidal systems. 

PACS numbers: 82.70.Dd, 82.70.Gs, 83.60.Rs 



I. INTRODUCTION 

Colloidal gels made of attractive particles at low vol- 
ume fraction raise both fundamental and industrial in- 
terest [J, [2l ■ For instance, among the huge variety of col- 
loidal systems, those made of fumed carbon black par- 
ticles suspended into an organic solvent or matrix are 
involved in many applications ranging from paints, coat- 
ings and printing to rubbers and tires 0] . Carbon black 
gels have also raised academic interest in the context of 
jamming 0, [H], shear-thickening Q, and dynamical 
heterogeneity as well as aging in attractive systems [1, Q . 
Besides their macroscopic rheological and aging proper- 
ties, recent research effort has focused on the effect of con- 
finement on such systems in search of new prospects in 
miniaturization and microfluidic applications. In partic- 
ular, it was shown that, when sheared within a confined 
geometry, various attractive particulate systems display a 
striking structuration into a striped pattern of log-rolling 
fiocs aligned along the vorticity direction Vorticity 
alignment was both inferred from light scattering experi- 
ments in colloid-polymer mixtures 11[ and Laponite clay 
suspensions and directly demonstrated by structure 
visualization in flocculated magnetic suspensions [l3l.[l^. 
non-Brownian carbon nanotubes p^ . attractive emul- 
sions p^ . and carbon black and alumina dispersions 
[l3 1 hinting to a somewhat general phenomenon. 

A mechanism for vorticity alignment was proposed 
based on an elastic instability of soft viscoelastic do- 
mains embedded into a less viscoelastic fluid: due to 
the streamline curvature inside the domains, an internal 
"hoop stress" develops that compresses the domains in 
the radial direction, leading to an elongation along the 
vorticity direction jisl This scenario is supported 

by the observation of small but noticeable negative first 
normal stress differences (Tsl - flTj . Although this effect, 



similar to a Weissenberg effect localized within viscoelas- 
tic domains, may indeed be involved in two-phase vis- 
coelastic fluids [1^, [1^ , it is less clear how compressible 
aggregates of attractive particles may play the role of vis- 
coelastic domains. Moreover, such a qualitative scenario 
does not predict the features of the striped pattern such 
as its wavelength and the width of the stripes. From the 
experimental point of view, the current literature suf- 
fers from a lack of detailed, quantitative characterization 
of the vorticity alignment phenomenon and only rather 
qualitative and fragmented data is available on the pat- 
tern static and dynamic features. In the particular case 
of carbon black gels, the vorticity-aligned structures were 
investigated only in rotational shear devices, such as the 
plate-plate and the cone-and-plate geometries, where ei- 
ther the shear rate or the gap width is spatially inhomo- 
geneous 0,[l3- Therefore, the aim of the present paper is 
to complement previous results on carbon black gels with 
a thorough analysis of vorticity-aligned patterns in vari- 
ous geometries, including simple plane shear, by varying 
all the control parameters (gel concentration, shear rate, 
and gap width). 

The present article is devoted to the steady-state char- 
acteristics of vorticity-aligned floes. In Sect. |TT] below, we 
first briefly explain the preparation protocol, recall the 
basic rheological characteristics of carbon black gels, and 
describe the translational and rotational setups used to 
shear the samples within gaps of width h = 40-400 fim. 
The experimental results are presented in Sect. lIIIl where 
we first qualitatively compare the patterns observed in 
three different geometries before turning to a quantita- 
tive image analysis of the influence of the various con- 
trol parameters. Finally our results, which should feed 
future modelling of vorticity alignment in attractive sys- 
tems, are discussed in Sect. [TVl in terms of compressible 
cylindrical floes. 
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II. MATERIALS AND METHODS 
A. Carbon black gels 

1. Sample preparation. 



clusters, leading to an increase of the effective particle 
volume fraction and therefore to enhanced viscous dissi- 
pation i.e. shear-thickening. From this shear-thickened 
gel state, vorticity-aligned structures are readily formed 
under steady shear at low shear rates 0, • 



Our gels consist of dispersions of carbon black (CB) 
particles (Cabot Vulcan XC72R) in a light mineral oil 
(Sigma, density 0.838, viscosity 20 mPa.s) at concen- 
trations ranging from 0.25 to 3 % w/w as described by 
Trappe et al. [8[ CB particles (density 1.8, typical diam- 
eter 200 nm, fractal dimension 2.2) are themselves small 
aggregates of permanently fused "primary" particles of 
diameter 20-40 nm. Our samples are prepared by dis- 
persing the CB powder into the oil through a vigorous 
manual agitation followed by a 1 hour sonication, which 
breaks any aggregates between CB particles. CB parti- 
cles interact through a short-range attractive potential, 
whose depth is estimated at about 30 kT for a 2 % w/w 
suspension Q. This attraction results in the formation 
of a weak colloidal gel that is easily disturbed by shear. 
Figured] shows the structure of ~30 /xm-thick samples as 
seen under a standard optical microscope. As previously 
shown by Trappe et al. [3, 01 , the gel is formed by a net- 
work of interconnected CB particles, which gets denser 
as the concentration is increased. 




FIG. 1: Optical microscopy images of carbon black gels at 
(a) 1 % w/w and (b) 2 % w/w. 



2. Rheological properties. 

The rheological measurements presented in this para- 
graph are performed at room temperature (~25°C) with 
a stress-controlled rheometer (Anton Paar, MCR 301) 
in a Plexiglas Couette cell of gap width 1 mm. In or- 
der to ensure a reproducible initial gel state, the sam- 
ple is presheared at a high shear rate (-1-1000 s~^ then 
-1000 s~^ for 20 s each) prior to any measurement. The 
gel is then left at rest for 100 s, which is sufficient to let 
the gel structure reform and reach a steady state with 
ne glig ible aging 

U,^- As discussed in previous studies 
[gI. Il7l| , preshear at high shear rates induces the breakup 
of locally dense clusters of CB particles into less dense 
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FIG. 2: Oscillatory shear rheology. Storage modulus, G' , 
(solid symbols) and loss modulus, G" , (open symbols) for CB 
gels at various concentrations: 1 % (o), 2 % (□), and 3 % w/w 
(A), (a) G' and G" as a function of the frequency / for a given 
stress amplitude of 0.2 Pa. (b) G' and G" as a function of the 
shear stress amplitude a for / = 1 Hz. Dashed lines indicate 
the points where G' = G", corresponding to the crossover 
from solidlike to fluidlike behaviour as a is increased. 



As shown in Fig. Ufa), the frequency-dependent vis- 
coelastic moduli of 2 % and 3 % w/w CB gels are typical 
of soft solids: the storage modulus G' is almost constant 
and a few times larger than the loss modulus G", which 
weakly increases with the frequency. Moreover, the stress 
sweeps of Fig. [2jb) allow one to roughly estimate the 
yield stresses of 2 % and 3 % w/w samples to 0.85 Pa 
and 2 Pa respectively. In the case of the very weak 1 % 
w/w sample, the crossover from solidlike to fluidlike can 
be distinguished around 0.2 Pa. Correspondingly, the 
frequency sweep of Fig. [DJa) performed at a stress am- 
plitude of 0.2 Pa shows that G" > G' for the 1 % w/w 
gel over the whole range of frequencies. Note that these 
estimates should be taken with care due to possible wall 
slip and time-dependent effects [l^l . Still these rheologi- 
cal data arc compatible with those found in the literature 
for similar concentrations [1, 0, [l^l . 
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B. Experimental shear cells 

1. Simple shear geometry. 

In a first series of experiments, we sfiall use a home- 
made simple shear cell that consists of two parallel glass 
plates. This setup is sketched in Fig. [3] The lower glass 
plate can be slowly translated using an endless screw. We 
use a stepper motor to control the rotation speed of this 
endless screw and therefore the applied shear rate 7. The 
gap width, h, is fixed using two thin aluminum spacers. 
The spacers (1 cm wide and 25 cm long) are glued over 
the whole length of the lower plate. A force is applied 
downwards on the upper plate through a metallic rod 
that is slighlty bent in order to force the contact between 
the spacers and the upper plate. By directly following the 
actual motion of the bottom plate using a video camera, 
we checked that this force is small enough so that the 
bottom plate moves at a constant speed wq with negligible 
stick-slip between the spacers and the upper plate over 
velocities down to 10 /xm.s^^. 



)L upper plate 
^ lower plate 
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shaker ^ 



I motor 



endless screw 



V sample 




FIG. 3: Sketch of the simple shear cell as seen from the side. 

As recalled above, the CB gel has to be presheared in 
order to get reproducible results and to reach a shcar- 
thickcncd state. To this aim, the metallic rod that keeps 
the upper plate against the spacers is attached to an 
electromagnetic shaker (LDS, V20I) that imposes a si- 
nusoidal displacement at a frequency of 10 Hz and with 
an amplitude of 3 mm. This corresponds to an oscillatory 
shear rate of amplitude ^ 1100 s~^. Before starting an 
actual experiment at low shear rate, such a large ampli- 
tude oscillatory shear is imposed for 50 s, after which the 
sample is homogeneously black over the whole surface 
of the glass plates. This indicates that any large scale 
heterogeneity in the gel structure has been broken up by 
shear. The sample is then left at rest for 10 s before the 
stepper motor controlling the motion of the lower plate 
is turned on. We checked that changing the preshear and 
rest durations does not affect our results. 

The gap width h — 173 ± 10 ^m was measured us- 
ing high-frequency ultrasonic echography and further 
checked by measuring the apparent surface covered by 



small oil droplets of known volume and squeezed between 
the two glass plates, h was found to be homogeneous 
over the whole surface of interest, a rectangle of typi- 
cal length 150 mm and width 45 mm. This home-made 
parallel-plate setup, hereafter referred to as translational 
shear cell, allows us to shear a very large surface of sam- 
ple {'-^IQ c\r?) with a homogeneous strain field over dis- 
tances as large as 100 mm, i.e. a total strain that can 
reach 7 = 60, 000 %, with shear rates ranging from 0.1 
to 3 s-^ 

In order to follow the evolution of the gel structure 
under shear, a CCD camera (Mikrotron MC1310) is fixed 
above the shear cell. This camera is used with a standard 
objective (Fujinon, focal length 16 mm) for wide-field 
imaging or together with a microscope objective (Leica, 
N Plan 10 x) to get an enlarged view of the samples (see, 
e.g.. Fig. The whole setup is lit by a LED backlight 
source (PHLOX, LEDW-BL-200x200) that provides a 
homogeneous white lighting over 200 mm by 200 mm. 



2. Rheo-optical setup. 

While the great advantage of our parallel plate setup 
is to provide a confined plane shear flow over more than 
10 cm, its obvious drawbacks are that (i) it does not al- 
low us to easily vary the gap width h and (ii) it does not 
give access to the shear stress a exerted onto the sample. 
Indeed tiny force measurements on a large plate dragged 
over large distances and with the constraint of keeping 
the gap width constant are very difficult if not impossi- 
ble. Therefore, in order to perform stress measurements, 
we turn to a rheo-optical setup similar to that used previ- 
ously by Osuji et al. @,0,[i3i although we shall not deal 
with the rheological measurements in the present paper. 

This rheo-optical setup, hereafter referred to as ro- 
tational shear cell, is sketched in Fig. |31 A rheometer 
(Anton-Paar, MCR 301) applies a shear rate 7 and mea- 
sures the corresponding stress response (j(t). This stress- 
imposed rheometer is used in a strain-controlled mode 
thanks to a feedback loop that allows one to impose a 
constant 7 within a time corresponding to a shear strain 
of less than 5 % for the range of shear rates investigated 
here. The CB sample is visualized from below through a 
mirror fixed under the transparent Plexiglas lower plate 
using the same camera and objectives as those of the 
translational setup. For the larger magnification, the mi- 
croscope objective is set between the lower plate and the 
mirror. 

In order to explore the influence of the geometry, we 
use either a Plexiglas plate (of diameter 60 mm) or a 
stainless steel cone (of angle 2°, diameter 50 mm, and 
truncation 209 ^m) as the upper rotating tool. Such a 
rotational setup does not suffer from any limitation in 
the accessible strain. The vertical position of the rotat- 
ing tool is controlled within 1 /zm. However, due to the 
tolerance in machining the bottom Plexiglas plate, the 
gap width h may vary locally by ±5 fim so that, for the 
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4: Sketch of the rotational shear ceU as seen from the 



plate-plate geometry, gap widths down to 40 fj.m only 
were considered. 

As recalled in the introduction, the plate-plate geom- 
etry has a homogeneous gap width but the stress and 
strain fields are strongly heterogeneous: the shear rate 
depends linearly on the radial position r as 7(r) = rflo/h 
where Qq is the rotation speed of the upper plate. In the 
following, wc shall denote by 7 the shear rate at the pe- 
riphery of the upper plate so that 7(r) = 'jr/R where R 
is the radius of the upper plate. On the other hand, the 
cone-and-plate geometry has quasi-homogeneous stress 
and strain fields but the gap width varies from h = 
209 fim at the cone apex to h = 875 /im at the cone 
periphery. 

The preshear protocol is as described above in 
Sect. Ill A 21 for rheological measurements. The only dif- 
ference is that, in the case of the plate-plate geometry, 
the rotation speed JIq (rather than the shear rate 7) is 
kept constant whatever the gap width in order to avoid 
expulsion of the sample for the largest gaps, fig is cho- 
sen so that 7 = 10, 000 s~^ for h = 100 fiui. Although 
this preshear protocol differs from the oscillatory proto- 
col used in the translational setup described above, it 
sets the gel in similar initial shear-thickened and repro- 
ducible states, allowing for a direct comparison between 
the various devices and geometries. 



III. EXPERIMENTAL RESULTS 

In this section results obtained with the setups intro- 
duced above arc described. Wc first give an account of 
qualitative observations made under simple shear and in 
both the plate-plate and conc-and-platc rotational ge- 
ometries. We then turn to a quantitative analysis focus- 
ing on the influence of the various control parameters, 
namely the sample concentration C, the gap width h, and 
the applied shear rate 7. In all cases, the data presented 
here arc taken in "steady-state," i.e. after the transient 
corresponding to shear-induced structuration and before 
any instability of the shear- induced structures |17| . A 



quantitative study of the initial transient and of long 
time stability will be the subject of future work. 



A. Qualitative observations 



1. Simple shear geometry. 




FIG. 5: Shear-induced structures under simple plane shear. 
The arrow denoted as vq shows the direction of shear. The 
dotted frame shows the typical area used for computing the 
wavelength. Inset: same shear- induced state at a larger 
magnification. These pictures were taken once the vorticity- 
aligned floes are fully developed for C = 1.5 % w/w, h = 
173 /im, and 7 = 1.34 s^^. 



Figure [5] shows the highly-ordered shear- induced state 
of the gel in the translational shear cell at two differ- 
ent magnifications. The dark line in the middle of the 
wide-field image corresponds to the shadow of the end- 
less screw. As already reported in a rotational setup by 
Osuji et al. the (opaque) CB particles arrange into 

long, dense floes aligned perpendicularly to the shearing 
direction. The oil is transparent and fills the gaps be- 
tween vorticity-aligned floes in the main picture. With 
the present parallel plate setup, straight floes with a di- 
ameter of about 100 ^m are observed that can be as 
long as 10 mm, which corresponds to an aspect ratio of 
100. The inset of Fig. [5] also shows that some large free 
clusters of CB particles remain in the interstitial space 
between the floes. Moreover, moviesf of the sample un- 
der shear at the largest magnification clearly show that 
the floes rotate as solid-bodies in a log-rolling motion. 
In such movies, it can be seen that in this shear-induced 
"steady-state," floes may break and recombine with other 
neighboring floes and that the oil flow in between the floes 
is probably much more complex than simple shear. 
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2. Rotational geometries. 

Figure [6] shows pictures of the gel after shear in the ro- 
tational plate-plate and cone-and-plate geometries. The 
qualitative features of the pattern of vorticity-aligned 
floes are in full agreement with previous observations 
In the plate-plate setup [see Fig. |6lja)], the wave- 
length seems to be constant in space but a large number 
of defects are seen due to the curved geometry. 




FIG. 6: Shear-induced structures in rotational geometries for 
C = 2 % w/w and 7 = 0.5 s"^ . (a) Plate-plate geometry with 
h — 250 fim. (b) Cone-and-plate geometry. The gray circle 
shows the limit between the conical and the plane region of 
the upper cone of truncation 209 fim. The white arrows point 
to conical floes. The white dots in (a) and (b) indicate the 
rotation axis. 

In the cone-and-plate geometry [sec Fig.lHl^b)], the two 
large white spots correspond to reflections of the inci- 
dent light on the metallic cone. Due to the rather large 
truncation of our cone, a circle was drawn to clearly dif- 
ferentiate between the inner plane region (for r < 6 mm) 
and the conical region (for r > 6 mm) of the tool. In 
this geometry, the shear-induced structures are wider at 
the cone periphery and thinner at the center. In some 
parts of the image, conical floes, whose apparent width 
increases with the distance from the center, can be clearly 
identified [see arrows in Fig. [S^b)]. Although this result 
may seem natural, this shows unambiguously that floes 
are not necessarily cylindrical and that their morphology 
is governed by the geometry of the shear cell. 
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FIG. 7: Spatio-temporal diagrams of (a) a horizontal section 
of Fig. [S] and (b) a circle of radius r = 20 mm in Fig. [HI a), 
(c) Advection speed Wstruct of the shear-induced structures 
normalized by the velocity vq of the bottom plate as a func- 
tion of the shear rate 7 in the translational setup for various 
concentrations (C — 1-2.5 % w/w) and h = 173 fj.m. (d) 
Rotation speed fistruct of the shear-induced structures nor- 
malized by the rotation speed Qq of the upper plate as a 
function of the radial position r in the plate-plate rotational 
setup for various shear rates (7 = 0.5-5 s~^), various gap 
widths [h = 100-400 fim), and C = 2 % w/w. The gray 
dashed lines correspond to 0.5. 



agrams consist of parallel black and white stripes that re- 
flect the motion of single structures at a constant speed. 
By measuring the slope of these stripes in either case, 
one can easily estimate the floe advection speed Wstruct 
[Fig- IZUc)] or their rotation speed fistruct [Fig. [Hd)]. In 
both cases, the speed of the structures is always half that 
of the upper plate whatever the shear rate, the position 
in the sample, the gap width, and the gel concentration. 
This result indicates a solid-body rotation of cylindrical 
floes between the two plates. As already noted above, 
this picture is directly conflrmed by moviesf recorded at 
a large magnification in the translational setup. 



B. Quantitative results from image analysis 

1. Advection speed of the vorticty-aligned floes. 

Figures[7|a) and (b) present spatio-temporal diagrams 
extracted from pictures taken in the translational shear 
cell and in the plate-plate rotational geometry during the 
steady-state regime. Such diagrams are obtained by plot- 
ting a horizontal section of the image in the translational 
case and a circle centered on the rotation axis in the ro- 
tational case as a function of time. The time origin t = Q 
is arbitrarily chosen in the steady-state regime. Both di- 



2. Wavelength and floe width extraction. 

The wavelength A of the shear-induced pattern is ex- 
tracted from the images shown above by Fourier trans- 
form. In the case of the translational shear cell, we first 
compute the 2D fast Fourier transform (FFT) I{kx,ky) 
of a given grayscale image I{x,y) over a small area of 
typical size 3 cm^ (see the dotted frame in Fig.O. This 
raw FFT is then filtered using a 2D convolution over 3 
points (7 resp.) in the ky [k^ resp.) direction in order to 
remove the noise [see Fig. [SJa)]. Finally, we average the 
2D FFT for -0.6 < ky < 0.6 mm"! and we look for the 
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maximum of (\I(k:c,ky)\) away from the origin and 600 

for the corresponding kx = kg to extract the wavelength 
A = 27r/A:o. Using the larger magnification (see inset of 
Fig. [5)), we also measure directly the apparent floe width 
d under simple shear. By selecting various parts of the 3 
sample, we checked that both A and d do not vary in '< 
space (data not shown). 
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FIG. 8: (a) Amplitude \I{kx, ky) \ of the 2D fast Fourier trans- 
form taken over the dotted frame shown in Fig. [5] and filtered 
as explained in the text, (b) Average of the 2D FFT shown 
in (a) over the fc^-axis for —0.6 ^ ky < 0.6 mm~^. The 
wavelength is readily extracted as A = 2n/ko, where fco is 
the wavenumber corresponding to the maximum away from 
the origin (see dashed lines). The FFT amplitude is coded 
in decibels with the reference taken at the maximum i.e. we 
plot 20 logjQ(|/(fci, fcy)|/ max(|I|). In (a), black corresponds 
to dB and white indicates levels smaller than -20 dB. 

In the case of the rotational setup, from images similar 
to Fig. [HI we extract the intensity /(r, 9) as a function of 
the angle at a fixed distance r from the rotation axis. 
Then we calculate the FFT of /(r, 9) over 9, whose max- 
imum yields the angular period that we multiply by the 
distance to the rotation axis to recover the spatial pe- 
riod A(r) of the pattern. As for the translational setup, 
the floe width d is extracted by directly measuring the 
apparent width of vorticity-aligned rolls from calibrated 
images obtained under large magnification. 

In all cases A and d are computed for 20 to 100 images 
recorded in the steady state regime. The values shown 
in the following figures are the average of these measure- 
ments and the error bars show the standard deviation 
of these measurements. As an example, the wavelength 
measured in the plate-plate rotational device for various 
values of r is presented in Fig. [5] These data show that 
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FIG. 9: Wavelength A of the shear-induced pattern as a func- 
tion of the radial distance r in the plate-plate rotational ge- 
ometry. C = 2 % w/w, h = 150 Atm, and 7 = 0.5 s"\ The 
dashed line shows the mean value of A over all values of r. 

A is spatially homogeneous so that what follows an ad- 
ditional spatial average over r = 10-25 mm will be used 
when analyzing the experiments performed in the plate- 
plate geometry. 

3. Influence of the shear rate. 

Figure [TU] presents the wavelength A and the hoc width 
d as a function of the applied shear rate 7 for various 
concentrations in the translational setup. Up to exper- 
imental uncertainty, A and d are independent of 7 over 
one decade in shear rate. This means that some shear 
strain is necessary to induce the structures but that the 
rate at which this strain is applied is not relevant to the 
steady-state regime. 

In the plate-plate rotational geometry, we recall that 
A does not depend on the radial distance r (see Fig. 
This also implies that A is independent of the local shear 
rate 7(7'). By systematically varying the plate rotation 
speed we found that A does not depend on the global 
shear rate 7, i.e. the shear rate at the periphery of the 
plate, either (data not shown). 

Moreover, there clearly exists a critical shear rate 7c 
above which structuration does not occur. We shall dis- 
cuss the measurements of 7c and the resulting "phase 
diagram" in more details in Sect. IIII B 61 below. 

^. Influence of the concentration. 

As seen in Fig. [TT] for the translational setup, both 
the wavelength and the hoc width do not show any clear 
dependence on the concentration C . The same result 
was obtained in the plate-plate rotational geometry (data 
not shown). Interestingly, this suggests that the mech- 
anism underlying the structure formation and selecting 
the pattern is insensitive to the gel concentration. It 
also clearly implies that the floe density has to increase 



7 



^ 400 





1 1.5 2 2.5 3 

C (% w/w) 



FIG. 10: Influence of the shear rate 7. (a) Wavelength A of 
the shear-induced pattern and (b) apparent floe width d vs 
7 in the translational shear cell {h — 173 /xm). The different 
symbols correspond to various concentrations: C = 0.5 % (o), 
1 % (□), 1.5 % (A), 2 % (.), 2.5 % (o), and 3 % w/w (■). 
The dashed hnes show the mean values of A and d over all 
values of 7 and C. 



FIG. 11: Influence of the gel concentration C. (a) Wavelength 
A of the shear-induced pattern and (b) apparent floe width d 
vs C in the translational shear cell {h = 173 /im). Each data 
point corresponds to an average over the various shear rates 
shown in Fig. 1101 The dashed lines show the mean values of 
A and d over all values of C. 



cides with the gap width ft,, as can also be checked in the 
with increasing concentration, which we shall discuss fur- translational setup in Figs [lUl HH and ^T[h) . 
ther in Sect. IIVI Above a concentration of 3 % w/w, CB 
gels were not observed to present any clear shear-induced 
structuration in our range of gaps and shear rates. 

6. Phase diagram for shear-induced structuration. 



5. Influence of the gap width. 

As shown in Fig. [T2T a) for a 2 % w/w gel in the plate- 
plate rotational geometry, A strongly depends on the gap 
width h. A linear fit yields A = 2.5h + 60 ^m. This linear 
scaling is confirmed in the inset of Fig. [T^] through ex- 
periments performed in the cone-and-plate geometry. In 
this geometry both the gap and the pattern wavelength 
are not uniform. However, when considering the local 
gap h{r) and the local wavelength A(r), one recovers a 
linear relationship A(r) = 2.0 h{r) + 50 /im, which is con- 
sistent with the previous ones. We conclude that the 
average linear law A ~ 2.25 ft. + 55 /im provides a good 
description of the data whatever the geometry, the gel 
concentration, and the applied shear rate. Still, due to 
experimental uncertainty, it remains unclear from these 
data whether the pattern wavelength tends to zero or not 
under extreme confinement ft — >■ 0. 

Finally, Fig. [T^ b) shows that the fioc width d also in- 
creases linearly with the gap width: d = 1.05 ft -I- 12 /im. 
Here, the y-intercept of 12 /im is believed to be insignifi- 
cant in view of the experimental uncertainty. Therefore, 
our results indicate that the apparent fioc width d coin- 



As already noted, shear-induced structuration docs not 
occur when the applied shear rate exceeds a characteristic 
shear rate jc- A simple way to measure jc is to take 
advantage of the non-uniformity of the shear rate in the 
rotational plate-plate setup, where the local shear rate 
7(r) spans the whole range from to the shear rate 7 
at the periphery of the plate. Indeed, if one imposes a 
large enough 7, one observes that vorticity-aligncd rolls 
develop only in the inner part of the sample while the 
outer part always remains homogeneous as seen in the 
inset of Fig. [T31 This allows for a precise measurement of 
the critical radius that separates structuration from no 
structuration. 7c is then defined as jc = jifc) ~ jfc/ R- 
We checked that different estimates of Tc extracted from 
images at various 7 at the periphery yield the same values 
of 7c to within 10 %. 

Figure [T3] shows 7c plotted in logarithmic scales as a 
function of the gap width ft for different concentrations. 
Here again, the data does not depend significantly on the 
concentration. The resulting curve, which separates the 
homogeneous regime from the structuration domain in 
the (7, ft) "phase diagram," is well fitted by a power-law 
behaviour 7c = Aft~" with a = 1.4 ± 0.1 
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FIG. 12: Influence of the gap width h. (a) Wavelength A of 
the shear-induced pattern vs h for C = 2 % w/w in the plate- 
plate rotational geometry. The solid line is A = 2.5 /i 4-60 ^im. 
Inset; A vs the radial distance r for C — 2 % w/w in the 
cone-and-plate geometry. For r > 6 mm, the solid line is 
2.0 h{r) + 50 /xm, where h{r) = rO is the local gap and 9 = 2° 
is the cone angle. For r < 6 mm, the gap width is constant and 
equal to h = 209 fj,m due to the cone truncation, (b) Apparent 
floe width d vs /i in the plate-plate rotational geometry for 
C = 2 % w/w (o). The black square (■) shows the floe width 
measured in the translational shear cell (h = 173 fim). The 
solid line is d = 1.05 h + 12 /im. 



IV. DISCUSSION AND CONCLUSION 

In summary, the main results of the present work are 
as follows. The steady-state features of shear-induced 
vorticity-aligned rolls are independent of the imposed 
shear rate 7 (provided 7 is smaller than a well-defined 
critical shear rate 7c) and of the concentration C of the 
CB particles. Whatever the shearing geometry, the only 
parameter that controls the shear-induced pattern is the 
gap width h and we found linear dependences d ^ h and 
A ~ 2.25 /i for the apparent floe width d and for the pat- 
tern wavelength A respectively. Finally, the critical shear 
rate jc is also independent of C and scales as h~^ '^. 

This last scaling may be explained in terms of a simple 
force balance between the viscous drag force acting on a 
floe and the attractive forces that keep individual parti- 
cles together. Such a force balance is classically used to 
estimate the typical cluster size as a function of the shear 
rate in flocculated suspensions under shear . Here, we 
consider floes generated from the inital gel network before 
structuration and elongation along the vorticity direction 
occurs. For confined CB gels, the typical size of such fiocs 
is given by the gap width h so that the viscous force can 
be estimated as i^visc ^ Vsjh'^ i where 77s is the oil viscos- 
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FIG. 13: Critical shear rate 7c, above which no structuration 
is observed, as a function of the gap width h and for various 
concentrations: 1 % (□), 1.5 % (A), 2 % (•), 2.5 % (o), and 
3 % w/w (■). The solid line is the best power-law fit of the 
fuU data set: > = A/i"" with A = 3650 and a = 1.39. 
Inset: determination of 7c in the rotational plate-plate setup 
for C = 2 % w/w, h = 300 fj,m, and 7 = 2.5 s~^ (see text). 
The white dot indicates the rotation axis and the white circle 
shows the critical radius beyond which no structuration is 
observed. 



ity. This viscous force will be able to tear the fioc apart, 
say in two pieces of size h/2, if it exceeds the attractive 
force Fattr ~ NgU/S, where 5 is the range of the attrac- 
tive potential, U its depth, and A''^ is the number of CB 
particles in a cross-section of the fioc. Introducing the 
floe fractal dimension df, one has TVs ~ h'^f~^ so that 
-Fattr ^ h'^f~^U/5. Therefore, structuration into floes of 
size h is predicted to be prevented by viscous drag for 



vise > ^attr; which corrcsponds to 7 > 7c with 

U 



7c 



(1) 



The critical shear rate is thus predicted to follow a 
power- law behaviour and, in the framework of Eq. ([l}, 
we deduce df = 3 — a = 1.6±0.1 from experimental 
observations. Previous works devoted to the equilibrium 
structure of different CB particles (including the ones 
used in the present study) suspended into various oils 
have reported fractal dimensions df = 1.8 ±0.1, hint- 
ing at a Diffusion Limited Cluster Aggregation (DLCA) 
mechanism [21I [22j . However, such an agreement be- 
tween the fractal dimension of the gel network at rest and 
that deduced from the experimental measurements of 7c 
should be mitigated since shear-induced fiocs are not ex- 
pected to keep the same fractal dimension as the intial gel 
state. Indeed, both experimental and theoretical studies 
of weakly aggregated dispersions have shown an increase 
of the fractal dimension from 1.8 to about 2.3 due to 
shear [23l - l25| . Therefore, the approximations leading to 
Eq. ([IJ may be too crude and a more subtle theoretical 
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(a) C = 0.25%1 

(b) (7 = 0.5^- 



is linked to the initial volume fraction of the gel at rest 

ip by 




FIG. 14: Shear-induced patterns under simple plane shear 
with h = 173 fim for different gel concentrations C. Left: 
large magnification. Right: low magnification. 



approach, such as that developed by Potanin,|26j would 
probably be better suited. In any case, it would be in- 
teresting to also test the influence of U and rjs, e.g., by 
adding a dispersant to lower the interaction potential. 

While the fact that the floe size is fixed by the gap 
width seems rather natural for confined log-rolling ag- 
gregates, a surprising result of the present study is that 
all the steady-state features of the vorticity-aligned rolls, 
namely A, d, and 7c, are independent of the gel concen- 
tration C. If one assumes a perfect pattern of cylindri- 
cal rolls of diameter h with a wavelength A, it is easily 
shown that the particle volume fraction within a roll (proU 



4A 4A poii ^ 
TTti irh pcB 



(2) 



where Poii and pcB denote the densities of the suspending 
oil and of the CB particles respectively and the last ap- 
proximation results from ip = PouC / [poiiC + /Ocb(1 — C)] 
with C <C 1. Since the wavelength A does not depend on 
C, Eq. ([2]) shows that the volume fraction inside vorticity- 
aligned rolls is simply proportional to C. Enlarged im- 
ages of the shear-induced structures at various concen- 
trations confirm qualitatively that the floes are very ten- 
uous at the lowest concentration [C = 0.25 % w/w, see 
Fig. im^a)] and get more and more compact as C is in- 
creased [see Fig. [HTb-c)]. 

Still, when taking a wider look at the shear-induced 
pattern (see right column of Fig. it clearly appears 
that the picture of an ideal striped pattern used in Eq. ([2]) 
is not very realistic outside the range C = 1-2 % w/w. 
Indeed, Fig. ll4f a-b) shows that lots of "holes" are present 
in the shear-induced pattern obtained at the lowest con- 
centrations (0.25-0.5 % w/w) while at the largest con- 
centrations (2.5-3 % w/w), very large agglomerates su- 
perimpose to the roll pattern [see Fig. [T^ f-g)]. This sug- 
gests that there exists a maximum volume fraction for 
the vorticity-aligned structures above which the rolls can 
no longer compress. However, in this case, it is unclear 
why the system does not lower the pattern wavelength 
to accommodate more rolls per unit area. Rather than 
a decrease in A, we observe a coexistence between very 
large clusters and the same roll pattern. 

To conclude, the present experimental data set paves 
the way for future modelling of shear-induced vorticity 
alignment in attractive systems. In particular, we be- 
lieve that theoretical works should concentrate on identi- 
fying the mechanism for structure formation and on the 
prediction of the pattern wavelength. Such a prediction 
may result not only from considerations on viscoelastic 
effects within the gel as already suggested in previous 
studies but also from a detailed analysis of the oil 

flow around the floes. Moreover, dynamical information 
on the pattern formation based on both image analysis 
and transient rhcology should provide even more experi- 
mental input to such theoretical attempts. 
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